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METHODS AND KITS FOR SCRgENlN£ UUCLEXC. I D_ DUPLE X STAB ILITY 
Introduction 

This invention was made in the course of research 
sponsored by the National Institutes of Health. The U.S. 
5 Government may have certain rights in this invention. 

Background of the Invention 

Mutagenic lesions in DNA frequently result from 
structural modifications of the heterocyclic bases (exocyclic 
adducts, free radical - induced modifications) and/or from 

10 complete removal of the base (abasic sites) . Further, 
cellular processing of lesion containing DNA can lead to 
mismatches, additions, deletions or base pair substitutions. 
Thus, by lesion, it is meant to include mismatches, additions 
and deletions. A wide range of lesion- induced thermodynamic 

15 effects have been observed. Typically, the free energy 
stabilizing the duplex is reduced significantly by inclusion 
of the lesion. Methods of determining lesion effects on 
duplex free energy are limited. Typically, the effect of a 
DNA modification on the energetics of duplex formation is 

20 measured by comparison of independently measured association 
constants for the modified and unmodified duplex or by 
comparing Tm values, which are commonly but erroneously 
believed to represent thermodynamic stability. Therefore, 
there is a need for a simple, reproducible and sensitive 

25 method for rapidly screening for duplex stability. 

The assays of the invention have two novel features 
which, when combined, provide a powerful and rapid method to 
assess the consequences upon duplex formation of perturbations 
to localized and/or global chemical features of nucleic acids. 

3 0 The first unique aspect is using two simultaneously competing 
equilibria to identify differences in equilibrium constants, 
between formation of two different nucleic acid duplexes. 
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Typically, the effect of a DNA modification on the energies 
of duplex formation have been measured by comparison of 
independently measured association constants for the modified 
DNA and the unmodified duplex requiring two separate 
5 experiments. A second novel feature of the present invention 
is that these assays require only one experiment. 

Furthermore, there is a large technical barrier for 
direct measurement of single duplex association events. In 
conventional titration experiments, a solution of one strand 

10 is added to a solution of its complement with formation of a 
duplex monitored by any of a variety of methods, including 
spectroscopic and calorimetric methods . To extract useful 
information from a conventional titration, the experiment must 
be devised such that a significant fraction of free titrant 

15 will be present throughout the titration. Satisfaction of 
this condition leads to the familiar sinusoidal shape of the 
titration curve. To satisfy this condition, typically the 
product of the initial titrate concentration, c, and the 
association constant, K, is in the range 10<cK<1000. Due to 

20 the high association constant for nucleic acid duplexes, the 
component concentration must be below the association 
constant, the components are likely to be too dilute to be 
detected by standard spectroscopic means. Having nucleic 
acids compete for duplex formation, as in the present method, 

25 creates a second equilibrium, referred to herein as a 
"competition" for duplex formation, that is measurable at 
essentially any concentration range. Thus, the concentrations 
can be tailored to virtually any method of detection. 
Further, the competition is measured directly from a single 

3 0 experiment rather than having to compare the results of two 
independently measured experiments. 

Another innovative aspect of this invention is the 
ability to discriminate between the two duplexes being formed. 
In one embodiment of the present invention fluorescence energy 

35 transfer is used to facilitate this discrimination. A common 
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spectroscopic method for monitoring duplex formation relies 
on the hyperchromicity of duplex formation. However, the 
extinction coefficients of duplexes of similar length is not 
a very sensitive reporter of the small differences in DNA 
5 content that are of most interest, as in the case of 
oligonucleotide duplexes with damage to only a single base. 
Even a technique such as circular dichroism is not 
sufficiently sensitive and also suffers from difficulties in 
interpretation of spectral variation which can be due to 
10 factors other than duplex formation. FET provides a unique, 
extremely sensitive, and essentially binary means of 
discrimination because only a duplex with both the donor and 
acceptor dyes will have the spectroscopic signature of the 
energy transfer. 

15 Competitive equilibrium assays of the present 

invention are more widely adaptable to a variety of nucleic 
acid systems than are assays that are based on changes in 
intrinsic spectral characteristics of the dyes. For example, 
the FET assay of the invention is dependent only on the 

2 0 presence of the two dyes and is limited only by the necessity 

of modifying DNA to bear the dyes and the fact that the 
distance between the dyes increases substantially when the 
initial FET duplex is disrupted. Any spectral changes that 
accompany the disruption of the FET duplex can be easily 
25 treated during data analysis and do not cause any significant 
complication for the FET assay. 

Thus, the present invention has a number of 
significant advantages over prior art techniques for 
determining duplex stability. 

3 0 Summary of the Invention 

An object of the present invention is to provide 
methods for screening for nucleic acid duplex stability by 
competitive equilibria. In these methods, a solution is first 
produced containing a known amount of an initial or reference 
35 nucleic acid duplex with a known stability. The initial 
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duplex is comprised of a first nucleic acid strand having a 
sequence, wholly or in part, homologous to a target strand and 
a second nucleic acid strand having a sequence, wholly or in 
part, complementary to the target strand. A series of 
5 additions of target strand are then made by titrating the 
solution with a second solution comprising a known 
concentration of the target nucleic acid strand. This target 
nucleic acid strand competes with the first nucleic acid 
strand for binding to a nucleic acid strand of the initial 

10 nucleic acid duplex. After each addition or titration, the 
solution is subjected to conditions which disrupt some or all 
of the nucleic acid duplexes and triplexes in the solution; 
subjected to conditions which promote duplex or triplex 
formation, and then monitored for any changes in the amount 

15 of initial nucleic acid duplex formed as a function of the 
amount of target nucleic acid strand added. This method can 
be used for extracting enthalpy data by controlling 
temperature during duplex or triplex formation and monitoring 
changes as a function of temperature so that a family of 

20 titration curves can be made and used to extract enthalpy 
(AH°) data. 

Another object of the present invention is to provide 
a for detecting a single nucleotide polymorphisms. In this 
embodiment of the invention, the initial nucleic acid duplex 

25 comprises a first and second nucleic acid strand, wherein the 
first or second strand of the duplex is designed to identify 
a single nucleotide polymorphism in a single- or double- 
stranded target nucleic acid sequence. In this method, the 
amount of the initial nucleic acid duplex in a solution is 

30 first determined. A fixed excess amount of a target nucleic 
acid strand is then added to the solution. The solution is 
then subjected to conditions which disrupt some or all 
duplexes or triplexes in the solution followed by conditions 
which promote duplex or triplex formation. The amount of 

35 initial duplex formed after addition of the target strand is 
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then measured. This measured amount, after addition of the 
target strand, is indicative of the target strand containing 
the single nucleotide polymorphism. 

Another object of the present invention is to provide 
5 methods for determining the concentration of a target nucleic 
acid strand which comprises adding a known volume and 
concentration of an initial nucleic acid duplex with a known 
stability to a known volume of a solution containing a target 
strand. Alternatively, a known volume of a solution of target 

10 strand can be added to a known volume of a solution containing 
a known concentration of an initial nucleic acid duplex with 
a known stability. The solution is then subjected to 
conditions which disrupt the initial nucleic acid duplex and 
any duplex between the target strand and a strand of the 

15 initial nucleic acid duplex followed by conditions which 
promote duplex formation. The relative change in the amount 
of initial duplex formed in the solution after addition of the 
target strand is used to determine concentration of the target 
strand . 

2 0 Another object of the present invention is to provide 

a method for assessing stability of various selected target 
strands. In this method competitive equilibrium assays are 
performed with the same initial nucleic acid duplex for each 
selected target strands. Changes in the amount of initial 

25 nucleic acid duplex formed as a function of the amount of the 
selected target nucleic acid strand added are compared for 
each target strand to ascertain differences in stability of 
duplexes or triplexes formed by the various target strands. 
In this embodiment of the invention it is not necessary to 

30 know the stability of the initial nucleic acid duplex. 

In a preferred embodiment of these methods of the 
invention, the first nucleic acid strand comprises a donor 
nucleic acid strand labeled with a donor of a FET pair and the 
second nucleic acid strand comprises an acceptor nucleic acid 

3 5 strand labeled with an acceptor of the FET pair and changes 
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in the amount of initial nucleic acid duplex in the titrated 
solution are monitored by measuring changes in FET donor or 
acceptor intensity. 

Yet another object of the present invention is to 
5 provide kits for screening for nucleic acid duplex stability 
and single nucleotide polymorphisms by competitive equilibria 
methods . 

Detailed Description of the Invention 

The present invention is a method and kit for 

10 screening the impact on nucleic acid duplex stability of 
alterations in base structures, e.g. due to carcinogen 
exposure or synthetic modification; the role of sequence 
context on such effects; and any other such alterations in 
Watson-Crick pairing such as mismatched base pairs or 

15 bulged/unpaired bases. Modifications of the phosphodiester 
backbone may also be detected. The small amounts of material 
required, speed of execution and applicability to a wide range 
of test strands make the assay useful for detailed 
thermodynamic characterizations and for screening 

20 applications. These features make the methods and kits of the 
present invention superior to calorimetry and other optical 
methods, with significant savings in the amount of test 
material necessary to perform a free energy determination. 
The range of effects on duplex stability of a particular 

25 defect or modification can be readily determined with the 
methods of the invention. With this information, the most 
interesting duplexes can be identified for further study. 
Thus, large investments in time and materials will be made 
only for systems of real interest. 

3 0 One of the key features of the invention is the 

flexibility in the nature of the test nucleic acid components 
that can be evaluated. In one experiment, two 13 mer DNA 
oligonucleotides, each bearing one of the FET f luorophores , 
referred to herein as donor (D) and acceptor (A) strands, 

3 5 form a duplex. The test or target strand competitor is a 
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third DNA oligonucleotide of the same length bearing a single 
damaged site at the central position. However , application 
of the method is general, encompassing virtually any variation 
in the nature of the nucleic acid duplex and test strand 
5 competitor. The donor (D) and acceptor (A) strands need not 
be the same length nor must the target, whether single strand 
or duplex- Further, as used herein, the terms "DNA", "nucleic 
acid" , "oligonucleotide" and "strand" are meant to include 
other variations as there is no requirement for any of the 

10 three nucleic acid strands to be DNA. Accordingly, the terms 
"DNA" , "nucleic acid" , "oligonucleotide" and "strand" are 
meant to include DNA, RNA, and analogues including those 
comprised, in whole or in part, of modified bases and/or 
modified backbones such as pept ido-nucleic acids (PNA) and 

15 other oligomers, incorporating modified phosphate and/or sugar 
moieties (e.g. PNAs , methyl phosphonates , phosphorothioates ) 
that maintain duplex- forming ability may be used in the method 
of the invention. Further, these terms are also inclusive of 
the vast number of non- Wat son- Crick nucleotide base variations 

2 0 that may be incorporated into any of the components, including 
intra strand crosslinks, abasic sites, naturally occurring or 
synthetic base variants, base mimetics and base adducts, 
including, for example, carcinogen- induced adducts. There is 
also no need to limit the system to three independent strands 

25 of the same length competing for formation of two possible 
duplexes . 

Most nucleic acid amplification techniques, such as 
polymerase chain reaction produce duplex target . The 
technique of the present invention can be used on such targets 
30 in one of two ways. 

In a first embodiment, a triple helix is formed 
between the target duplex and one of the strands of the 
reference duplex. The sequence requirements for triple helix 
formation are well known. Typically, triple helices involve 
35 stretches of pyrimidines on one strand of a Watson-Crick 
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duplex, a complementary stretch of purines on the other strand 
of the Watson-Crick duplex, and a third strand comprised of 
stretches of either complementary purines or pyrimidines which 
resides in the major groove of the duplex. Details of the 
5 sequence requirements and tolerated substitutions are well 
known and widely described in the art . The target duplex need 
not be disrupted and all calculation equations provided in the 
Examples are applicable to this embodiment without 
modification . 

10 The second embodiment is more complex and applicable 

only when sequences do not meet the requirements for triple 
helix formation. In this embodiment, the target duplex must 
be disrupted, i.e. melted, so that both the donor- labeled and 
acceptor-labeled oligonucleotide can bind to the complementary 

15 strand of the target duplex. The addition of two equilibria 
(the formation of the target duplex and the interaction of the 
donor strand with one of the target duplex strands) makes this 
formalism inappropriate for deriving quantitative data. 
However, due to the law of mass action, the equilibrium 

20 distribution of complexes will still depend on the relative 
values of the equilibrium association constants and the 
various concentrations. Therefore, the FET observable will 
also depend on these values. As a consequence, qualitative 
information on the relative stability of complexes formed by 

25 the target duplex's component strands and various sets of DNA 
duplex probes can be obtained. 

The method is also useful for structures that might 
include large bulging/unpaired regions, competing internal 
loops and/or hairpins or other deviations from simple duplex 

30 formation. The above-mentioned variants may occur in 

combination, thereby increasing the number of potential 
targets of study. 

The only limitations are that the FET donor and 
acceptor be within resonance distance in the initial duplex 

35 and that formation of the competing complex prevents energy 
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transfer by displacement of either donor or acceptor. Any FET 
donor and acceptor pair can be used and such dyes are well 
known in the art and commercially available. The fluorescent 
dyes may be at opposite ends of the duplex {-5' and -5' or 3 ' 
5 and 3'), the same end of the duplex (-5* and -3 1 ) , or with one 
or both fluorophores in the interior of the strand (s) . The 
fluorophores may be linked after oligonucleotide synthesis or, 
when the phosphoramidites are available, incorporated during 
synthesis . 

10 Further, FET, monitored either by fluorescence of the 

acceptor or by quenching of the donor, is not the only usable 
means of monitoring the amount of the reference or donor- 
acceptor duplex. Any method by which the amount of reference 
or initial duplex can be monitored as a function of the amount 

15 of target can be used. Eximer fluorescence or other optical 
means can be used. In fact, nucleic acid strands of the 
initial duplex can be labeled with any pair of species with 
properties or characteristics dependent upon proximity. 
Examples include, but are not limited to, fluorescent dyes, 

20 antibody-antigen pairs, enzyme- inhibitor pairs and enzyme- 
coenzyme pairs. Further, if the assay is carried out on a 
surface, surface plasmon resonance (SPR) spectroscopy may be 
employed with the label being a chromophore at the wavelength 
used in the SPR measurement . 

25 The high throughput nature of this assay, in 

comparison to the more time and material intensive techniques 
generally used for thermodynamic analysis, makes the assay 
applicable to various problems in biotechnology and 
pharmaceutical research. For example, this assay can be used 

30 to evaluate nucleotide mimetics as drugs such as the anti-HIV 
drugs ddC and AZT; to evaluate the effects of 
carcinogen/chemical exposure on the stability of DNA and DNA- 
RNA hybrid duplexes; and for screening of various parameters 
for hybridization studies (e.g. temperature, buffer, 

35 sequences) . Screening of non-natural nucleic acid analogs as 
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antisense or antisense agents can also be addressed by this 
method. The assay of the invention can be a companion 
technique to help improve existing hybridization techniques. 
The method of the invention can also be used for screening, 
5 in solution, for the presence of single nucleotide 
polymorphisms (SNPs or "snips") which are used in 
pharmacogenetic research targeted at identifying the genetic 
basis of disease and genetic diagnosis of the potential for 
such disease in individuals. The latter aspects have 

10 particular importance in the biotechnology industry. Many 
companies are currently involved in developing and marketing 
hybridization assays for a wide variety of research and 
development efforts. Virtually all of the assays currently 
in use rely on immobilization of at least one participant in 

15 the hybridization reaction. Significantly, the immobilization 
introduces a host of complications, including non-specific 
interaction of any/all of the components with the immobilizing 
platform; possible distortion of the biochemically important 
equilibrium due to immobilization; the possibility that the 

20 chemical linkage of the immobilization can partially occlude 
the necessary interactions with the non- immobilized 
components; and the necessity of additional steps in the 
protocol for the immobilization itself prior to running the 
binding experiments. The method of the invention, being 

25 entirely in solution, eliminates the complications caused by 
immobilization. Further, because the method uses titration, 
control experiments with standardized DNA can be run 
frequently or in parallel with test compounds to eliminate 
spurious results. Such standardized DNA is a component of a 

3 0 kit for carrying out the method of the invention. 

The competing equilibria which are the basis of this 
assay provide a greatly enhanced method for detecting 
differences in stability between two nearly identical 
duplexes. Studying the association of two strands forming one 

3 5 duplex and comparison of the association of two other strands 
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in a separate experiment, as is done in current methods, 
requires two experiments with the inherent compounding of 
experimental error. The present invention advantageously 
requires only one. 
5 The calculations used for data analysis are derived 

from the general equations for three component, two-equilibria 
systems as taught by Linn and Riggs (J. Mol . Biol. (1972), 72, 
671-90) . 

Two equilibrium association constants are defined 
10 below, with subscripts f and t indicating free and total 
concentrations, respectively, and AD, AX, D, A, and X 
represent the donor/acceptor complex, the competitor/acceptor 
complex, donor strand, and acceptor strand, and the competitor 
or "test" strand, respectively: 



K \J2C \ \JX 



iOT _ (w t >j)H:>j>j) 



15 (Equation 1) 



K - lJD 



M7"|Dj.-l^J)([Aj,-[^J>j) 



(Equation 2) 



These equations are for the "test" or "target" strand forming 
a complex with the acceptor strand. In this and the models 
that follow, the opposite competition, with X binding to D, 
20 can also be described by simple substitution of the terms. 

These basic expressions for the equilibrium constants 
can be combined and rearranged to the following equation: 



WO 00/37686 



- 12 - 



PCT/US99/30751 



(Equation 3) 

The loss of energy transfer is monitored as the AD 
complex is disrupted by the formation of AX. The 
complementary measurement of emission of the acceptor 
5 subsequent to energy transfer also may be used. Which 
approach works better will depend on the photophysical 
properties of the fluorescent dyes. As discussed, the value, 
9 is the fraction of the initially observed fluorescence 
energy transfer at each point in the titration and is related 
10 to the relative concentrations of the donor/acceptor complex 
([AD]) to [total donor strand] ( [D] t ) ([AD] = [D] t at the 
beginning of the titration) : 

(Equation 4) 

with relative fluorescence readings for the dilution-corrected 
15 relative fluorescence at each point n, I n/ fluorescence of the 

fully formed FET duplex, i^, and the fluorescence of the donor 

strand alone, I D . 

In equation 5, the dilution corrections are added 

explicitly, where V 0 is the initial volume of the D-containing 
20 solution, V A , the volume of the added A-containing solution, 

and v 2 , the volume of the ith aliquot of the X-containing 

solution . 
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V + V + v 



I I 



V 



e 



V +v 



I -I 



V 



{Equation 5) 



Substituting Equation 3 into Equation 4, assuming 
that [X] t >> [AX] , and rearranging yields Equation 6. 



e 




(Equation 6) 

5 A program (written in the Microsoft VBA language) that 
computes a function theta ( [A] t , [D] t , [X] t , K AD , K^) to 
calculate the isotherms (G vs. X t or logX t ) using equation 6 
is depicted in Example 5. The value of 0 is found by 
iteration to satisfy the equation RHS (equation 6) - 0 
10 0, where RHS means right-hand side. The entire 

experimental isotherm can be fit using equation 6 to find a 
value for the desired parameter . This is, however, not 
necessary as shown in the following section. 



15 strand X at which 0 = 0.5, or exactly half of the 

acceptor/donor duplex, AD, has been disrupted. The value 
of X 0 5 can be interpolated from a plot of 0 versus log [X] t . 
When 0 = 0.5, a simple relation between the desired 
equilibrium constant K AX , the measured value X 0 s , and the 

20 known values [D] t and K AD results in: 



•0 . 5 



is defined as the concentration of competing 
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K = 



* K \d] 



2X 



(Equation 7) 



Application of the well known relation between AG° and K 
yields Equation 8: 



AG * =-ROi 



IX 



(Equation 8) 

5 The value of K M will have been previously determined by 
independent methods such as differential scanning calorimetry 
and UV-absorbance melting. Alternate representations of 
equation 8 can be used to evaluate the free energy changes 
associated with the formation of the duplexes and the defects. 

D 

AG ° -AG ° -XT In ■ — 

2X 



10 (Equation 9) 



D 

AG ° =~RT InK -RT In - 

2X 



(Equation 10) 
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AG." -AG : =-RT In- 



2X 



{Equation 11) 

Alternatively, the impact of a difference between two 
oligonucleotides (X 1 & X 2 ) on duplex stability, AAG°, can be 
evaluated by titration (in separate experiments) of the two 
5 competitors against the same reference AD duplex at the same 
AD concentration (Equation 12). 



A AG" , =AG.; -AG/ =-KT ii 



~X~ 



(Equation 12) 

In Equation 12, the value for K M is canceled. Therefore, one 
can evaluate the free energy impact of a single defect, or 
10 multiple defects, without thorough thermodynamic 

characterization of the AD duplex. 

The assumption of [X] t >> [AX] in equation 6 can be 
relaxed thereby leading to 

g . K(i-e) 

X K 1 *K(i-e) 



(Equation 13) 

15 An unknown parameter [AX] must be evaluated as part of the 
calculation of 6. This is accomplished by iteration over [AX] 
with the restriction that [A] t - [AX] + [AD] + [A] f . A second 
root finding problem is executed with the objective being 
finding a value of [AX] that satisfies the equation 
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[a] -\jk]-bd f- ,=0 

' tf, (1-0) 

(Equation 14) 

A second program is also shown in Example 5 that calculates 
theta( [A] t# [D] t/ [X] c , K^, K M ) without the restriction that 
[X] t >> [AX] . While a useful reduction of equation 13 
5 (similar to equations 7 and 8) cannot eliminate the need for 
iterative solution, values of X c at G = 0.5 are readily 
calculated. Values of X t so calculated are compared to 
experimental values and K M adjusted to produce the 
experimental value of X t . 

10 Comparison of isotherms calculated using equations 

6 and 13 reveal that the error introduced by the assumption 
of negligible [AX] is small and only significant when K M ~K M . 
For K^/K^ - 1, the error in free energy is about 0.4 
kcal/mole; for K M /K^ = 0.1, 0.06 kcal/mole; and for K^/K^ = 

15 0.01, 0.005 kcal/mole. Therefore, the assumption is in most 
cases reasonable. Those cases where it is not totally without 
adverse consequence, that is where K M -K M( are readily 
predictable, can be addressed easily by application of the 
full equation 13 . 

2 0 As is clear from the equations 7 and 8, the K M and 

AG° M values that are measured are relative to the value of 
K^. As a practical matter, titration of X into the solution 
to a concentration of 1000 [D] t is convenient. An X 0 5 value of 
1000 [D] c corresponds to a AAG° value of about 4.5 kcal/mole 
25 (AAG° - -RTln ( 1/2000) . This is a rather large range and 
should accommodate most single base defects . The range can 
easily be extended by performing additional titrations using 
less stable AD duplexes. The stability of the AD duplexes can 
be modulated by inclusion of modified bases and/or mismatches. 

3 0 A series of AD duplexes can therefore be designed to cover 
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essentially any range of AG° values, in intervals of 3 to 4 
kcal/mole . 

Analysis of the parallel titrations can be 
accomplished independently or collectively to determine the 
5 free energy values using the analysis described above. 

Equation 6 can be rearranged and used to determine 
the concentration of a target sequence (X) , when the values 
of and K^x are known. In this example, a known volume and 
concentration of AD duplex is added to a known volume of an 
10 X containing solution. Alternatively, a known volume of X 
containing solution is added to a known volume of AD duplex 
containing solution of known concentration. Thus [A] t = [D] t is 
known. The concentration of X, [X] t , can be calculated from 
the relative change in fluorescence, 9, using the formula 

[x] t =^{[A] t a-er -e}/e 



15 (Equation 15) 

When the temperature is controlled during the 
annealing process of a titration, additional information can 
be obtained. The cooling process can be performed in steps 
20 so that values of 9 are collected as a function of 
temperature. Data at each temperature are used to produce a 
family of titration curves. Each curve is analyzed 

independently and values for K M are determined as a function 
of temperature. The van' t Hoff equation, 

AH " =-R(d3rK /d(l/T) ) 



25 (Equation 16) 
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can be used to extract enthalpy ( AH° ) data. The thermodynamic 
description, at a given temperature is completed by using 
AG° = -RTlnK and AS°=(AH° - AG°)/T. 

Direct detection of formation of DNA duplexes by 
5 titration is extremely difficult because of the very low 
concentrations required for monitoring the equilibrium which 
are well below the operating range of traditional in-solution 
methods. A competition assay is usable over a wide range of 
instrumentally accessible concentrations. The method of the 

10 invention provides a more reliable measurement of nucleic acid 
complex stability over a very wide range of free energy values 
because the titration depends on the difference in stability 
between the initial donor/acceptor-containing duplex and the 
resulting competitor- containing duplex and not on their 

15 absolute free energy values. Therefore, the range of 

accessible free energy values can be tuned by choice of the 
initial donor/acceptor-containing duplex. Relative free 
energies are usually the desired experimental result and the 
method of the invention provides them directly. The lower 

20 detection limits of the fluorophores define the maximum 
difference in free energy that can be detected by the method. 
The use of fluorophore detection provides great sensitivity. 
The emission spectrum of the donor strand at 10 pM 
concentration has been visualized reliably using a photon- 

2 5 counting fluorometer. 

Free energies calculated from the assay of the 
present invention have been demonstrated to be in agreement 
with those measured by extensive thermodynamic studies on 
individual duplexes. In these experiments, two titrations 

3 0 were performed at the same D t concentration, for two starting 

Watson-Crick FET duplexes, designated A»T and T«A, which 
differ only by the central base pair, out of the 13 pairs in 
the duplex. Competition on each duplex is from nearly the 
same single strand as present in the FET duplexes, except this 
3 5 single strand is unlabeled and has a tetrahydrof uranyl abasic 
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"lesion" site (F) at the central base pair. Free energy 
values measured by this method compare quite favorably to 
those measured by extensive differential scanning calorimetry 
and UV absorbance melting experiments on these 13-mer duplexes 
5 containing a single tetrahydrof uranyl abasic site (F) in the 
central position. Specifically, using the FET assay, a value 
of -14.5 ± 0.1 kcal/mole for formation of the F-T duplex was 
determined, compared with a value of -15.1 ± 0.6 kcal/mole 
determined using DSC/UV melts. Similarly, using the FET assay 

10 a value of -16.2 ± 0.1 kcal/mole was determined for formation 
of the F*A duplex compared with a value of -16.0 ± 0.4 
kcal/mole by DSC/UV melts. These results correspond to AAG° 
values of -1.7 ± 0.2 kcal/mole from FET and -0.9 ± 1.0 
kcal/mole from DSC/UV melting studies for substitution of an 

15 A residue for a T residue opposite the abasic site. 

Because measurement depends on relative 
concentration, [D] J [X] , the method of the present invention 
can be used in any concentration regime. The appropriate 
concentration range is determined by the sensitivity of the 

20 detection system. Accordingly, [D] t is selected by one of 
skill to optimize the detection method. Due to practical 
limitations on the volume of titrant that can be added to the 
titrate, there are some practical limitations to the range of 
free energy values that can be covered by a single titration. 

25 However, this range is quite large. A convenient limit of the 
ratio [D] t /[X] 05 is about 0.001. This corresponds to a factor 
of about 0.002 in association constant or 3.7 kcal/mol in free 
energy. Most defects for which reliable free energy data are 
available fall into this range. 

3 0 However, since the measured free energy values depend 

on the ratio of and { AAG° ) , the appropriate choice of 

donor- acceptor duplex must be made for each titration. 
Because the complementarity of the component strand of the 
donor-acceptor reference duplex need be only sufficient to 

35 form the duplex and achieve moderate thermal stability, the 
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assay may be tuned over a wide range of free energy values by 
introducing mismatches in the donor- acceptor reference duplex. 
By judicious choices of donor-acceptor duplexes, a series of 
three donor acceptor duplexes can cover a range of 11 
5 kcal/mole in free energy. If titrations are done in parallel, 
use of such a family of duplexes relieves one of estimating 
the magnitude of prior to performing the titration. In 

this embodiment of the method of the present invention, 
several donor/acceptor complexes are formed simultaneously. 

10 Each different donor or acceptor oligonucleotide differs 
slightly, while maintaining complementarity of the acceptor 
with the target (X) or employing an analogous system where X 
binds to the donor, to produce donor/acceptor pairs of 
differing stability. By appropriate selection of donor and 

15 acceptor dyes, it is possible to conduct multiple simultaneous 
titrations of X into a single solution. Spectroscopic 
discrimination of the various donor acceptor pairs provides 
multiple free energy determinations simultaneously. 

The assay of the present invention can be adapted by 

20 immobilization to a variety of inert solid supports using 
technologies established for standard hybridization studies. 
In this embodiment, one of the strands of the reference duplex 
(either donor (D) or acceptor (A)) can be attached to the 
surface- The reference duplex can then be formed on that 

25 surface. Exposure to target will release the unattached 
strand, thereby producing the signal. Because concentration 
cannot be defined at a surface in the same way as in solution, 
comparison to a standard of known stability is required for 
quantitative results. However, immobilization facilitates the 

30 miniaturization and adoption of this assay to high throughput 
screening while providing the benefit of using FET and the 
competing equilibria to increase the sensitivity of 
measurement of differences in duplex stability. 
Immobilization may also allow for the construction of an 

35 array of donor/acceptor duplexes of differing stability. Such 



WO 00/37686 PCT7US99/30751 

- 21 - 

an array assures appropriate selection of the initial complex. 
Further, because ranges of the free energy differences 
measurable using the various initial complexes overlap, 
multiple complementary measurements can be made 
5 simultaneously. 

Simultaneous titrations provide a number of 
additional advantages to this assay. With appropriately 
designed donor-acceptor reference complexes, the range of 
accessible free energy is multiplied by the number of duplexes 

10 titrated simultaneously. Thus, employing three simultaneous 
donor-acceptor duplex titrations means that the effective 
range of a single titration experiment becomes 9-12 kcal/tnol, 
without expending any extra test strand. The enhanced range 
also means that a favorable outcome is likely in a single 

15 experiment, rather than having to explore various single donor 
acceptor duplexes to find one which has a free energy less 
than 3-4 kcal/mol higher than the test duplex. Second, this 
application provides a degree of multiplexing that reduces the 
time involved in each assay, and thereby enhances the ability 

20 to perform high throughput screening of nucleic * acid 
variations. Thus, the simultaneous titrations are truly 
simultaneous rather than merely in parallel, meaning that all 
of the information is gathered using a single cuvette. The 
detailed theory supporting the simultaneous titration 

25 experiment is described below. 

Single Acceptor/ Multiple Donor or Multiple Acceptor /Single 
Donor Method 

In principle, any number of donor-acceptor pairs 
could be included with selective excitation of the donors and 

3 0 a single acceptor. The limitation is imposed by the necessity 
of finding several donors with non-overlapping excitation 
spectra and sufficient Stokes shifts such that their emission 
spectra each overlap sufficiently the excitation spectrum of 
the acceptor. This method allows a series of AD complexes 

35 with varying K M to be used simultaneously. 
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As an example, the case wherein there are 3 donors 



(D 1# D 2 , D 3 ) and a single acceptor dye (A) is described; 
however, the number of donors is not limited to 3 . The donors 
are discriminated by the excitation wavelengths, Aj , A 2 , A 3 . 
5 The initial concentrations are [DJ t = [D 2 ] t = [D 3 ] t = [A] c /3. 
The Donor and acceptor labeled strand can form any or all of 
the complexes AD X/ AD 2 , AD 3 , and AX. The equilibrium 

constants for the various complexes that form are enumerated 
below . 



K 




(hj.-h A J)(N.4 D . A J>. A J-p. A J>J) 
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10 (Equation 17A) 
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(Equation 18) 



Analogously, three theta values, which differ by the 
15 excitation wavelength, can be defined. 
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(Equation 19C) 
and thus 



A D 



D 



I 3 -J 3 

0 3 " 3 P 3 



JD he Id 



5 (Equation 20A) 




(Equation 20B) 
(Equation 20C) 

Combining equations 20 with equations 17 and 18 and knowledge 
of the three values, allows for fitting for the value of 
10 . In principle, not all values need be known. In this 
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case, the increased number of unknowns complicates the 
analysis significantly. 

An alternate model, in which AD lf AD 2 and AD 3 are in 
equilibrium with XD X/ XD 2 and XD 3 can be derived analogously. 
5 Fitting here is more complex as the number of unknowns is 
larger - including 

K 

XD 2 

K 

XD 3 

A treatment in which a single donor and multiple 
acceptors can be used when fluorescence energy transfer can 
be observed directly (when acceptor emission is observable) 
10 and the acceptor emission spectra do not overlap 
significantly. Derivation of appropriate equations for this 
case is straightforward and analogous to those derived for the 
above case . 



Multiple Acceptor /Multiple Donor Methods 

15 Several donor acceptor/complexes can be monitored in 

solution simultaneously. Each donor must have a unique 
excitation wavelength and each acceptor an absorbance spectrum 
corresponding to the emission of its donor. If emission of 
the acceptor is to be monitored, the emission spectra of the 

20 acceptors must be unique. Some correction for overlap can be 
made; however, such necessity complicates the data analysis 
significantly. The number of donor acceptor pairs is, in 
principle, unlimited. 

Multiple donor acceptor pairs are designated A 1 D 1/ 

25 A 2 D 2 , A 3 D 3 , etc. An example is provided using three AD pairs, 
but any number is possible. Equations accounting for the 
multiple simultaneous equilibria are described below. 
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(Equation 21A 
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(Equation 21B) 
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(Equation 21C) 

|AX| _ 

^ 1 1 " W, [A j f ~ (W . "[A X\ -[A ~[A Xj)([A [ -[A £> j-[A D j-[A D ]-[A Xj) 



A X 



(Equation 22) 

5 It is assumed that each X sequence will bind to A 1 , A 2 , and A 3 
with equal affinity, since the three acceptor oligonucleotides 
have identical sequences. Therefore, [XAJ = [XA 2 3 = [XA 3 ] and 



K =K =K 



1 3* 



(Equation 23) 

Similar reasoning leads to the assertion that [DiAj] = [D^] 
10 = [D^] and so forth for the other donor sequences. Defining 
expressions for 6 and substituting the following expressions 
for the equilibrium constants it is found that: 
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(Equation 24A) 



[Equation 24B) 



(Equation 24C) 
Therefore , 



5 (Equation 25A) 



(Equation 25B] 
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13 = AD=AD=AD 



(Equation 25C) 

Further assume that [d } ] = [z> 2 ] = [z) 3 ] ' wnicil is determined 
the experimental setup. 
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(Equation 27) 



5 Equations for 



and 



K 



* 2' 



K 



assume similar forms and as noted above are assumed equal to 



K 
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This system of equations (Equations 26 & 27) must be 
solved by iterative methods to find a value of K M which 
satisfies, at each X t , the constraints imposed by the known 
concentrations of the AD complexes derived from the measured 
5 6 values. 

A fundamentally different strategy for simultaneous 
monitoring of multiple donor/acceptor pairs is to use the dyes 
as acceptor and donor, but on different duplexes. Again, 
discrimination is made optically. Here, an example is 
10 provided using three donor/acceptor complexes, but the method 
is not limited in the number of complexes that can be 
employed . 

The nomenclature is altered slightly from that used 
above. Here, each fluorescent dye is designated as D, with 

15 superscript A indicating that dye D is acting as an acceptor 
and superscript D indicating that dye D is acting as a donor. 
Dyes are attached to oligonucleotides such that three FET- 
capable duplexes can form: D A 2 D D 1 , D A 3 D D 2 , and D A 4 D D 3 . Dye D 1 acts 
only as a donor and D 4 only as an acceptor; however, dyes D 2 

20 and D 3 act as acceptor and donor, but on different duplexes. 
Again, the oligonucleotides are designed so that D A 2 D D 1I D A 3 D D 2 , 
and D A 4 D D 3 vary in stability systematically and so that the X 
strand can form duplexes with the acceptor bearing strands, 
namely D A 2 X, D A 3 X and D A 4 X . 

25 As in the cases described above, equilibrium 

constants can be derived relating the concentrations of the 
various solution components. 
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WO 00/37686 



PCT/US99/30751 



- 29 - 



-\d*d»\-\d«d a 



(Equation 2 8B) 



-\D 4 A D° 



(Equation 28C) 

Because the acceptor strands are not identical (the dyes 
differ although the oligonucleotides to which they are 
5 attached are identical) , there are three additional 
equilibrium constants to define. 
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(Equation 29A) 
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(Equation 29B) 
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(Equation 29C) 

.0 If the dyes do not perturb the equilibria significantly, it 
is reasonable to assume that °* x ^ D ^ x ^d 4 a x ^ The 
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assumption that the acceptor dyes are equally perturbing or 
non-perturbing is reasonable and testable. 
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5 (Equation 30C) 

O^ 1 6^ 

Note that in the expressions for n and " , the measured 

quantities {shown in parentheses) contain contributions from 

the fluorescence of the acceptor of the previous {as in the 

assigned index) donor acceptor pair. This is assumed to be 
10 independent of the formation of the complex. This assumption 

is testable, should it not be valid appropriate corrections 

can be applied. 

As in the above analyses, the concentrations of the 

donor/acceptor complexes can be determined by measurement of 
15 the G values and knowledge of the total concentrations of the 

donor strands. 
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(Equation 31C) 

Again substituting values for the 0 terms and assuming that 
5 = [A D ] = \^°\ ' W ^ 1 ^ C ^ L ""~ s determined by the experimental 

setup, a system of equations is presented that can be solved 
as a function of [X] t to find values for 
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(Equation 33) 
Expressions for 



K 



5 and 



are derived similarly and, as described above, their values 
are assumed to be identical to 

K A 

a 2 * 



The advantage of this method of multiple simultaneous 
titration, where some of the dyes act as both donor and 
0 acceptor (although on different oligonucleotides) , over the 
method described above, in which all donors and acceptors are 
unique, is that the probability of identification of a set 
of fluorescent dyes with the necessary photophysical 
properties is enhanced. 
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While there are restrictions on the properties of a 
set of dyes usable in simultaneous titrations, there are 
likely to be many sets of usable dyes. One such set of dyes 
is described in the table below. Each of them is available 
5 from Molecular Probes, Inc. of Eugene, Oregon. 
Donor Acceptor Fluorescent Dye 



D? 



Alexa 350 
Alexa 43 0 
Alexa 4 88 



D A 

10 4 Alexa 594 

When one of the reference duplex strands is attached 
to a surface, only the other strand need be labeled. In this 
embodiment, the measurement is not FET, but rather 
fluorescence. A washing step must be included in this 

15 embodiment to remove released label. However, this 

modification may disturb the equilibria rendering this 
embodiment of the assay more qualitative than quantitative. 

Multiple measurements can be conducted simultaneously 
on a surface. A series of different strands can be attached 

20 to a surface so as to identify a particular oligonucleotide 
with its location on the surface. This kind of spatial 
distribution of different oligonucleotides is well known. The 
fluorophore can be attached post - synthetically or as a 
phosphoramidite ; either method is compatible with methods for 

25 producing an oligonucleotide array on a surface. A single 
oligonucleotide with sufficient complementarity to form 
reference duplexes with each of the immobilized strands can 
be used to form a series of reference duplexes . The surface 
is then exposed to a single target that simultaneously 

30 equilibrates with all of the reference duplexes. A series of 
standards can be routinely incorporated into this assay. 
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The method of the present invention is of particular 
use in single nucleotide polymorphism (SNP) screening. The 
theory behind the application of this FET assay to SNP 
screening is based upon two segments of DNA , the first being 
5 a "wild type" sequence (does not contain a SNP) , and the 
second being a variant sequence (a SNP) that, for example, may 
be a marker for disease or disease tendency. Standardized 
methods abound for amplifying a small genetic sample (e.g. 
from a few microliters of blood) , and the result would be one 

10 amplified strand, which would become the unlabeled competitor 
(X) strand in our FET assay. When amplification of the target 
is performed the quantity of the target can be determined by 
use of labeled primers. Either fluorescent ly labeled, so as 
not to interfered with the subsequent FET measurement, or by 

15 any of the many well known methods for labeling amplification 
products . 

Since SNP screening is merely a detection of a SNP, 
and not a quantitation of energetic impact per se, FET-based 
SNP screening does not even require a full titration of target 

20 (X) into the donor- acceptor duplex. Merely adding a fixed 
excess amount of the X strand is sufficient, in effect making 
a two-point titration wherein the fluorescence is measured 
before and after X addition. The amount of X strand needed 
should be in the range of 10- to 100- fold excess over donor 

25 (D) . In these experiments, when the sequences of X and D 
match, X will efficiently displace D from the donor acceptor 
duplex, causing an efficient, and probably complete, reduction 
in 6. When X and D do not match, 6 will remain high. In 
theory, either experiment would be sufficient to prove the 

30 presence or absence of wild-type sequence. However, it is 
generally preferable to guarantee that a positive and negative 
signal will be achieved in clinical assays, so it would be 
preferable to run two experiments in parallel, monitoring 
donor-acceptor duplexes of wild- type and SNP sequences, one 

35 of which should have high 0 and one low 6. Additionally, a 
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heterozygous sample would score as an efficient competitor for 
both donor-acceptor duplexes. 

The assay is further applicable to situations where 
there may be multiple SNP sites within a single gene, with 
5 each variation occurring at low frequency. In this embodiment, 
a much longer sequence can be screened similarly, still in a 
single experiment. In such a case, the donor-acceptor duplex 
has the wild-type sequence of the full length of the region 
of interest. The amplified target strand, X, is judged, from 

10 a single 10- to 100-fold addition, to either compete or not 
for the donor-acceptor duplex. If FET is reduced, and X is 
a good competitor, then further screening could be done using 
additional assays more specific in nature to the sequences in 
question. Even if another specific screening method is 

15 preferred, the global nature of the FET assay is useful as a 
first screen, eliminating the cost and expense of screening 
by other more detailed methods for every single patient being 
tested. If, among all SNPs identified within a gene, the 
frequency runs as high as 20% for having any one variation, 

20 the single screen using a very long sequence would still 
eliminate the need for extensive testing on 80% of the 
samples. This would represent a large savings in material, 
and would greatly improve the throughput of testing 
facilities . 

25 The SNP assay is not limited by the length of the 

strands (e.g. number of bases long) that are being screened. 
The actual length limitation is a theoretical barrier 
attributable to a kinetic situation whereby very long uniquely 
complementary sequences may not ever form a complete duplex. 

30 This limitation is similar to the fact that long sequences 
such as complete bacterial genomes cannot be completely 
reannealed after thermal disruption. The barrier for this 
assay is estimated to be in the range of several hundreds of 
bases, such that sequences of interest in the range of up to 

35 100 bases can be easily accommodated. Clearly, longer 
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sequences, however, are considered within the scope of this 
assay . 

The relevance of a strand length of about 100 bases, 
however, is that regions bearing even one SNP may be assessed 
5 by a single assay, using this method. As SNP identity becomes 
linked to disease potential in humans, such screening becomes 
an extremely valuable technique for rapid diagnosis of 
individual patients, both for diagnosis of ailment, infection, 
drug sensitivity, drug resistance, etc., and early detection 

10 of conditions that can lead to timely application of 
preventative treatments. Regions of genes already identified 
as having relevant SNPs are generally in the range of 100 
bases or less, meaning that the theoretical length limitation 
will not become a factor in practice. 

15 The key feature of this method that eliminates length 

of the oligonucleotides as a limitation, unlike in other 
assays under commercial development, is that this method 
depends on competition between two different equilibria 
characterized by two different equilibrium constants. Thus 

20 the actual measurement is the ratio of the two equilibrium 
constants, and thus the difference in free energies for 
formation of these alternative duplexes. The difference in 
free energy caused by a defect, such as a single base pair 
mismatch (i.e. not a canonical Watson-Crick pairing), is the 

25 same whether the defect is within a very long duplex or a 
short one, so long as the bulk of the duplex generally remains 
stable. Thus, the magnitude of AG° is not a consideration, 
since the fluorescence assay measures the difference in two 
AG° values (AAG°) directly. Mismatches of the four canonical 

30 bases will generally exhibit AAG° values of about 1 or 2 
kcal/mol, but this 1 to 2 kcal/mol is the actual magnitude of 
the measurement itself. Thus, whether the absolute values are 
20 kcal/mol of duplex, as in a 10-15 base pair 
oligonucleotide, or 600 kcal/mol of duplex, in something 
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perhaps as long as a gene, the assay is silent to these 
magnitudes, because it is reporting AAG° directly. 

Fluorescent ly labeled strands used in the present 
invention may be provided in a kit form. A researcher 
5 interested in adapting their test DNA lesion could do so 
simply by synthesizing a small amount of test DNA within a 
sequence predefined by the kit DNA . A kit would include at 
least one FET- labeled starting duplex and the appropriate 
buffer along with simple instructions for performing the 

10 experiment and analyzing the data. The kit may also include 
a standardized competitor strand to ensure reproducibility and 
to facilitate comparisons over time. As an example, the 
strands of a kit could be designed with specific targets in 
mind. A specific kit could be designed, for example, to 

15 screen for a single mutation in a gene. 

The following nonlimiting examples are provided to 
further illustrate the present invention. 
EXAMPLES 

Example 1: Purification and dye-labeling of 5 1 -amino- linker 

2 0 oligonucleotides 

Standard phosphoramidite DNA synthesis with an amino- 
linker phosphoramidite as the last (5') residue is performed. 
The intact synthesis column is dried under vacuum. The column 
is then cracked open and the glass support beads are 

25 transferred to a screw top plastic bottle. Aqueous NH 4 OH (1 
ml) from the freezer is added to each tube and the DNA is 
deprotected for three days at room temperature for standard 
amidites . Tubes are then cooled in the freezer and the NH 4 OH 
is pipetted off. The supports are then washed with 2 x 200 

30 /il water or a mixture of EtOH : CH 3 CN : H 2 0 (3:1:1); the samples 
and washes are combined; and then dried in a Speed-Vac. The 
dried samples are then dissolved in H 2 0 at low temperature 
(<40°C) and floating non-soluble materials are removed. At 
this step, the sample may be purified by reverse phase, using 

3 5 high performance liquid chromatography (HPLC ) and a PRP- 
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column, equilibrated with 50 mM ammonium bicarbonate. Elution 
is performed by a linear (5-50%) gradient of acetonitrile in 
50 mM NH 4 HC0 3 . Following f reeze-drying, the fractions 
containing the purified tritylated N-modified oligonucleotide 
5 are heated to 95°C for 5 minutes to remove the MMT- (trityl) 
group and subjected to ethanol precipitation in the presence 
of sodium ions, as follows: the volume is adjusted to 300 to 
400 /il water; NaCl (100 jil , 2M) or sodium buffer is then added 
along with 950 fil of EtOH; samples are then placed in the 

10 freezer for at least 30 minutes, preferably one hour. 
Following freezing, the sample is centrifuged for 12 minutes 
at 14,000 rpm. The resulting pellet is dried to remove any 
residual ethanol. Trityl is then removed by addition of 200 
fil of 80% acetic acid for one hour at room temperature and 

15 evaporating the liquid in a Speed -Vac for 2 to 3 hours until 
a glassy residue is observed- This step is critical to ensure 
elimination of any residual free amines that might interfere 
with the labeling reaction. For the unlabeled 

oligonucleotides, the trityl group is removed by addition of 

20 80% acetic acid and incubation for one hour at room 
temperature, followed by freeze drying. Both the purity of 
the final product and the success of detritylat ion are 
monitored by analytical reverse phase HPLC. If required, 
additional purification of the detritylated oligonucleotide 

25 is performed by reverse phase HPLC. The purified DNA is 
subjected to ethanol precipitation/Na+ exchange to rid the 
sample of any excess reactive amines from the HPLC buffer and 
labeled as follows: DNA (30-40 OD-260) is dissolved in 270 fil 
of H.O in an O-ring tube. NaHCO 3 (30 /il, 1 M) at pH 8 . 3 is 

3 0 then added . One mg of succinyl ester form of the dye is then 
added for each 2 0 OD of DNA. This can be weighed as a dry 
reagent into a glass vial, dissolved in 80 /xl/mg of fresh 
DMSO, and added into the plastic tube of DNA solution. The 
glass vial is then rinsed with 20 of additional DMSO and 

35 added to the plastic tube. Alternatively, a 100 fil aliquot 
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of dye is added. The dye and DNA are then allowed to react 
at 3 7°C or higher at least overnight. 

The DNA is then isolated from any unreacted dye with 
a PD-10 Sephadex G-25 column (Pharmacia) . The column is 
5 equilibrated by rinsing with at least 25 ml of H 2 0 . The DNA 
sample is diluted with H 2 0 to a final volume of 1000 fxl and 
loaded onto the column. The DNA is then washed using 1 . 6 ml 
H 2 0. DNA is eluted in 600 /il fractions with H 2 0. Generally 
six fractions are sufficient. 

10 The fractions are dried to at least X A the volume and 

adjusted to 300 jxl total volume with H 2 0 . Each fraction is 
then independently subjected to ethanol precipitation. 
Generally only the first three fractions will have appreciable 
DNA. The free dye stays mostly dissolved in the ethanol. 

15 The DNA containing precipitates are combined and the 

labeled DNA purified using ion exchange HPLC (Mono-Q column) 
with Buffer A being 50 mM Tris HCl with 15% acetonitrile and 
Buffer B being Buffer A with 1 M NaCl . The gradient profile 
can be tailored but in general increases from 0 to 80% buffer 

20 B over 25-35 minutes. Monitoring absorbance at both 260 nm 
and the absorbance maximum of the attached fluorophore can 
help to define the labeled and unlabeled DNA peaks. 

The labeled DNA is further purified by HPLC using an 
ion exchange column (e.g. Mono Q, Pharmacia) equilibrated with 

25 50 mM Tris HCl and 15% acetonitrile (Buffer A) and a linear 
gradient (0-100%) of Buffer B (i.e., Buffer A containing 1 M 
NaCl) . Absorbances at 26 0 nm and the wavelength corresponding 
to the maximum absorbance of the fluorophore are use to 
monitor and define labeled and unlabeled pools of 

30 oligonucleotides. 

Alternatively, a desalting column can be used in 
place of the ethanol precipitation steps. 

5' -labeled oligonucleotides forming duplexes have 
been compared with the parent unlabeled duplexes, revealing 

35 no alterations in their thermodynamic stability in the 
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presence of the probes. Moreover, the ability to form Watson- 
Crick duplexes in a stoichiometric ratio has been further 
confirmed by HPLC analysis of the annealed mixtures, in 
comparison to the free oligonucleotide strands. Additional 
5 evidence reveals that both labeled and unlabeled 
oligonucleotides may be recovered upon completion of a FET 
assay by repurif ication, with no indication of adulteration. 
Time-dependent studies of the integrity of the labeled 
oligonucleotides reveal no sign of aggregation or degradation 

10 within 1 year from sample preparation. The labeled 

oligonucleotides may be stored as stock solutions in water and 
working buffer at or below -20°C for at least six months. 
Preferably, the samples should be stored as a lyophilized 
powder, for periods exceeding six months. 

15 Example 2: Determination of labeled DNA concentration 

Determination of the concentration of the labeled DNA 
strands in stock solutions has been performed using an average 
extinction coefficient of 1.1 x 10 5 M^cm" 1 at 25°C. The 
intrinsic DNA absorbance at 260 nm has been demonstrated to 

2 0 not be significantly altered by the presence of the conjugated 

dye . 

Example 3 : Formation of the FET duplex 

A working stock solution of each of the labeled DNA 
of complementary sequence in the range of 2 to 3 micromole DNA 
25 strand is prepared. Fluorescence detection is performed as 
follows: after each aliquot (of any titrant) the cuvette is 
heated to about 90 °C, using an external heat block, and cooled 
to 25°C via the intrinsic cooling of the jacketed cuvette 
holder in the fluorometer. The cuvette must be tightly 

3 0 stoppered to minimize evaporation during heating. Wavelengths 

must be tailored for the dyes used. For example, for the 
fluorophore pair Oregon Green 514 and Rhodamine -Red-X , the 
emission spectrum is collected over 510-650 nm with excitation 
at 508 nm, with scanning at 100 nm per minute. The time 
35 drive is set to collect, using the kinetics mode, 30 or 60 
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seconds (at 0.1 second per reading) using 508 nm excitation 
and 528 nm emission. These data points are then averaged, 
resulting in a precise relative fluorescence intensity for 
each reading, and an associated standard deviation for that 
5 averaged value. The fluorescence of the buffer alone is read 
to establish a blank for the instrument response. 

The fluorescence of the "free" donor strand is then 
determined. A sample of 10 nM donor strand is prepared in 250 
fil total volume of buffer and fluorescence is measured. An 
10 aliquot of the acceptor strand from the working stock 
sufficient to achieve a final concentration of 100 nM is then 
added to form the FET duplex and the fluorescence is 
determined . 

Example 4: Titration of the competing strand 

15 A working stock of the competing strand is prepared 

at a concentration appropriate to the expected 
ability/inability to compete for duplex formation with the 
formed donor/acceptor pair. In general, the concentration is 
one order of magnitude higher than the donor and acceptor 

20 working stocks for each 1 kcal/mol of free energy difference 
(AAG) expected for the competing strand. In the event the 
free energy differences are completely unknown, a number of 
concentrations can be prepared, covering a wide range of 
concentrations, i.e., 2 /iM, 20 /iM, 200 fxM , etc. solutions in 

25 buffer. Titration is started with the most dilute solution 
and more concentrated solutions are used as necessary. 

For the first aliquot, the competing strand is added 
to a final concentration of about ^ of the concentration of 
the donor/acceptor concentration. The fluorescence is then 

30 determined. If there is a significant change in the 
fluorescence, titration is continued with the working stock. 
If there is no change, titration is continued with a higher 
working stock concentration. Additional aliquots are added 
and the fluorescence determined until the fluorescence 
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intensity recovers to at least % of the intensity measured for 
the free donor strand. 

Example 5: Automation of competitive titration experiment 

The FET data acquisition protocol was automated on 
5 an AVIV Model ATF-105 Automatic Titrating Fluorescence 
Spectrophotometer. The customized software developed for the 
FET assay on this particular instrument improved the overall 
accuracy, precision, and data throughput compared to 
conventional manual titration experiments. There are several 

10 key features of the automated FET assay including programmed 
titration of acceptor and competitor or target strands and the 
ability to conduct successive heating/cooling cycles of the 
sample solution in the cuvette. Selection of the upper 
temperature limit is dictated by two important considerations, 

15 namely that the parent and test duplexes are dissociated into 
single strands and the covalently attached fluorophores are 
stable at the desired upper temperature. Extensive control 
studies conducted on the stability of the Oregon Green 514 
(OG) and Rhodamine-Red-X (RdRX) labeled strands indicate that 

20 these fluorophores retain their integrity during successive 
heating/cooling cycles over the temperature range of 0 - 75°. 
The superimposition of UV melting/ cooling curves reveled that 
the f luorescently labeled strands in the parent duplex are not 
labile below 75°C. The overall reproducibility is compromised 

25 when heating the identical duplex to temperatures above 75 °C, 
as noted in the family of non- superimposable UV 
melting/cooling curves. It is incumbent upon the analyst to 
judiciously evaluate and select the practical upper 
temperature limit for a particular duplex and set of 

3 0 fluorophores. 

A typical experiment is initiated by placing a 
cuvette containing a 100 nM solution of the Oregon Green 514 
labeled donor strand in the sample compartment, heating the 
cuvette to 75°C, maintaining the temperature at 75°C for three 

35 minutes, and cooling the sample to 20°C over an equilibration 
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period of five minutes. During the cooling cycle, the 
instrument is operated in the kinetic mode and the 
fluorescence emission intensity is recorded at 521 nm (i.e. 
excitation wavelength = 508 nm) at five second intervals. 
5 After the fluorescence intensity plateaus indicating that 
equilibrium is achieved, a reference spectrum of the donor 
strand is recorded over the wavelength range of 510-650 nm. 
Activation of the first syringe drive dispenses fixed 
microliter aliquots of the Rhodamine -Red-X labeled acceptor 

10 strand in the sample cuvette to form the FET-active donor- 
acceptor reference duplex. Upon addition of each aliquot of 
acceptor, the sample is subjected to a heating/cooling cycle 
identical to that above to ensure that the duplex anneals 
properly. The donor strand (D) d (GCGTACACATGCG) -OG (SEQ ID 

15 NO:l) is titrated with its complementary acceptor strand (A) 
d ( CGCATGTGTACGC) -RdRX (SEQ ID NO : 2 ) to form the FET-active 
reference donor-acceptor duplex. The diluted corrected 
fluorescence intensities are cast in the form of a Job Plot 
to confirm that the stoichiometry of the single strands in the 

20 reference duplex is 1:1 (i.e. mole fraction = 0.5). Although 
the forward titration demonstrates that the single strands 
associate to form a competent duplex, the experimental 
protocol may be simplified by loading the pre- formed reference 
duplex into the cuvette prior to conducting the automated 

25 competition experiment. Elimination of the forward titration 
increases the sample throughput by reducing overall 
experimental time in the FET assay by approximately 50 
percent . 

Having formed the FET active reference duplex in the 
30 initial part of the experiment, the second syringe drive is 
activated to dispense fixed aliquots of the unlabeled 
competing strand ( X ) d ( CGCATGFGTACGC ) (SEQ ID NO : 3 ) . The sample 
solution containing the three strands is subjected to 
heating/cooling cycles after each addition in the titration 
35 experiment to facilitate annealing of the donor-acceptor and 
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acceptor- target (AX) duplexes. A sufficient excess of the 
target strand X (in this case approximately 100 fold) is 
titrated into the cuvette to ensure that at least half of the 
acceptor has been displaced from the reference duplex. The 
5 dilution corrected relative fluorescence (6) is plotted as a 
function of the concentrations of acceptor (A) and target (X) 
strands. The concentration of competing strand (X 0 . 5 ) at which 
exactly half of the acceptor/donor duplex (AD) is disrupted, 
that is 6 = 0.5, is interpolated from this plot. Substituting 
10 the value for X D b into the simplified relation: 

K M = (K AD -D t /2«X 0 . 5 ) 

yields a value for the acceptor/target association constant 
(K M ) . Application of the thermodynamic relation between AG° 
and facilitates calculation of the free energy: 

15 AG°=-R-T«lnK AX 

In this example, values of X os = 1.1 x 10" 5 for 
d ( CGCATGFGTACGC ) ( SEQ ID NO : 3 ) and D t = 4.22 x 10" 8 for 
d ( GCGTACACATGCG) -OG (SEQ ID NO : 1 ) , coupled with a 1^=9.0 x 
10 14 determined independently for the reference duplex, results 

2 0 in a x 10 12 for the formation of the AX duplex. 

Substitution into the relevant relation yields a value of AG° 
= 16.4 kcal/mol for the AX duplex that compares favorably with 
the value determined previously (i.e. AG° = 16.0 kcal/mol) 
employing a combination of calorimetric and spectroscopic 

2 5 technique. 

Example 6 : Equation Programs 
Equation 6 Program 
'declare globals ' 

Dim At, Dt, Xt, Kad, Kax As Double 
30 Function theta(zl, z2 , z3 , z4 , z5) As Double 

At = zl 

Dt = z2 

Kad = z3 

Kax = z4 
35 Xt = z5 
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theta = ModRegFal (0#, 1#) 
End Function 

Private Function f unci (A) As Double 

fund = At * (1 - A) / {(1 + Xt * Kax) / Kad + Dt * (1 - A) 
5 - A 

End Function 

Private Function ModRegFal (XI , X2 As Double) As Double 
'Modified Regula Falsi' 

1 (adapted from Conte & de Boor, 1980) ' 
10 'Finds root of function fund, if bracketed 1 

Const xtol = 0.000000000001 
Const ftol = 0.000000000000001 
Const ntol = 100 

Dim SignFl, n, PrvsF3 As Integer 
15 Dim Fl, F2, F3 , X3 As Double 

Fl = fund (XI) 
F2 = funcl(X2) 

'test whether root is bracketed 1 
If Sgn(Fl) * F2 > 0 Then 
20 Debug. Print "Xl= ", XI, "X2 = X2 

Debug. Print "fund (XI )= ", Fl , "funcl(X2)= ", F2 

End 
End If 

X3 = XI 
25 F3 = Fl 



1 BEGIN ITERATION' 
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'TEST FOR CONVERGENCE ' 

' Is interval small enough? 1 
If Abs(Xl - X2) <- xtol Then 
5 ModRegFal = X3 

Exit Function 
End If 

' Is F3 small enough? ' 
If Abs(F3) <= ftol Then 
10 ModRegFal = X3 

Exit Function 
End If 

1 GET NEW GUESS BY LINEAR INTERPOLATION 1 
X3 = (Fl * X2 - F2 * XI) / (Fl - F2 ) 
15 PrvsF3 = Sgn(F3) 
F3 = funcl(X3) 

•CHANGE TO NEW INTERVAL ' 
If Sgn(Fl) * F3 >= 0 Then 
XI = X3 
20 Fl = F3 

If F3 * PrvsF3 >= 0 Then F2 

Else 

X2 = X3 
F2 = F3 

25 If F3 * PrvsF3 >= 0 Then Fl 

End If 



= F2 / 2# 



= Fl / 2# 



Next n 

1 END ITERATION' 
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Debug. Print "Xl = " , XI, M X2= ", X2 , M X3= M , X3 

Debug. Print "fund (XI )= Fl , "funcl(X2)= ", F2 , "funcl(X3)= 
" , F3 

Debug . Print ntol, " iterations without convergence" 
5 End Function 



Equation 13 Program 
'declare globals • 

Dim At As Double, Dt As Double, Xt As Double 

Dim XA As Double, Kad As Double, Kax As Double 
10 Dim AX As Double 

Function theta(zl#, z2#, z3#, z4#, z5#) As Double 

At = zl 

Dt = z2 

Kad = z3 
15 Kax = z4 

Xt = z5 

dummy = ModRegFal2 { 0# , zl) 'switch to At 
theta = ThetaOO 

Debug. Print "AX/Xt= "; AX / Xt ; " Theta = "; theta 
20 End Function 

Private Function ThetaO ( ) As Double 
ThetaO = ModRegFall (0#, 1#) 
End Function 

Private Function f unci {A As Double) As Double 
25 fund = At * (1 - A) / ( (1 + (Xt - AX) * Kax) / Kad + Dt * (1 
- A) ) - A 
End Function 

Private Function func2 (B As Double) As Double 
Dim TO As Double 
30 AX = B 

TO - ThetaO ( ) 
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func2 = At - AX - TO / (Kad * ( 1 - TO ) ) - TO * Dt 
End Function 



Private Function ModRegFall (XI As Double, X2 As Double) As 
Double 

5 'Modified Regula Falsi 1 

1 (adapted from Conte & de Boor, 1980) 1 

'Finds root of function fund, if bracketed 1 

Const xtol = 0.000000000001 
Const ftol = 0.000000000000001 
10 Const ntol = 100 

Dim SignFl As Integer, n As Integer, PrvsF3 As Integer 
Dim Fl As Double, F2 As Double, F3 As Double, X3 As Double 

Fl = fund (XI) 
F2 = fund (X2) 

15 'test whether root is bracketed' 

If Sgn(Fl) * F2 > 0 Then 

Debug. Print "root not bracketed" 

Debug. Print "Xl= ", XI, "X2 = " , X2 

Debug. Print "fund (XI ) = Fl, " funcl(X2)= " , F2 
2 0 End 
End If 

X3 = XI 
F3 = Fl 

'BEGIN ITERATION 1 
25 For n = 1 To ntol 
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1 TEST FOR CONVERGENCE 1 

' Is interval small enough? ' 
If Abs<Xl - X2) <= xtol Then 
ModRegFall = X3 
5 Exit Function 

End If 

'Is F3 small enough?' 
If Abs(F3) <= ftol Then 
ModRegFall = X3 
10 Exit Function 

End If 

'GET NEW GUESS BY LINEAR INTERPOLATION ' 
X3 = (Fl * X2 - F2 * XI) / (Fl - F2 ) 
PrvsF3 = Sgn(F3) 
15 F3 = fund (X3) 

' CHANGE TO NEW INTERVAL ' 
If Sgn(Fl) * F3 >= 0 Then 

XI = X3 

Fl = F3 

20 If F3 * PrvsF3 >= 0 Then F2 

Else 

X2 = X3 
F2 = F3 

If F3 * PrvsF3 > = 0 Then Fl 
25 End If 

Next n 

' END ITERATION 1 
Debug. Print "XI = " , XI, "X2= ", X2 , "X3 = " , X3 



= F2 / 2# 



- Fl / 2# 
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Debug. Print "fund (XI) = " , Fl, "funcl(X2)= F2 , "funcl(X3)= 
" , F3 

Debug -Print ntol , " iterations without convergence" 
End Function 

5 Private Function ModRegFal2 (Yl As Double, Y2 As Double) As 
Double 

'Modified Regula Falsi' 

' (adapted from Conte & de Boor, 1980) 1 

'Finds root of function f unc2 , if bracketed ' 

10 Const xtol = 0.000000000001 
Const ftol = 0.0000000001 
Const ntol = 100 

Dim SignFl As Integer, n As Integer, PrvsF3 As Integer 
Dim Fl As Double, F2 As Double, F3 As Double, Y3 As Double 

15 Fl = func2 (Yl) 
F2 = func2 (Y2) 

'test whether root is bracketed' 
If Sgn(Fl) * F2 > 0 Then 

Debug. Print "Yl= " , Yl, "Y2= ", Y2 
20 Debug. Print "func2(Yl)= ", Fl , "func2(Y2)= ", F2 

End 
End If 

Y3 = Yl 
F3 = Fl 



25 1 BEGIN ITERATION' 
For n = 1 To ntol 
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' TEST FOR CONVERGENCE ' 



■Is interval small enough? 1 
If Abs(Yl - Y2) <= xtol Then 
ModRegFal2 = Y3 
5 Exit Function 

End If 



'Is F3 small enough?' 
If Abs(F3) <= ftol Then 
ModRegFal2 = Y3 
10 Exit Function 

End If 



' GET NEW GUESS BY LINEAR INTERPOLATION 
Y3 = (Fl * Y2 - F2 * Yl ) / (Fl - F2 ) 
PrvsF3 = Sgn(F3) 
15 F3 = func2 (Y3) 



1 CHANGE TO NEW INTERVAL' 
If Sgn(Fl) * F3 >= 0 Then 

Yl = Y3 

Fl = F3 

20 If F3 * PrvsF3 >= 0 Then F2 = F2 / 2# 

Else 

Y2 = Y3 
F2 - F3 

If F3 * PrvsF3 >= 0 Then Fl = Fl / 2# 
25 End If 



Next n 

•END ITERATION 



Debug. Print "Yl= " , Yl, "Y2 = ", Y2 , "Y3= " , Y3 
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Debug. Print "func2(Yl)= ", Fl, "func2(Y2)= ", F2 , "func2(Y3) 
" , F3 

Debug. Print ntol, " iterations without convergence" 
End Function 



